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Abstract. We study the long-time behavior of infinite-energy solutions to the incom- 
pressible Navier-Stokes equations in a two-dimensional exterior domain, with no-slip 
boundary conditions. The initial data we consider are finite-energy perturbations of a 
smooth vortex with small circulation at infinity, but are otherwise arbitrarily large. Us- 
ing a logarithmic energy estimate and some interpolation arguments, we prove that the 
solution approaches a self-similar Oseen vortex as f — ^ oo. This result was obtained in 

■ collaboration with Y. Maekawa (Kobe University). 

Comportement asymptotique en temps des solutions de V equation 
Q . de Navier-Stokes dans un domaine exterieur du plan 

D 

Q- - - 

' , Resume. Nous etudions le comportement asymptotique en temps des solutions de I'equa- 

■ tion de Navier-Stokes incompressible dans un domaine exterieur du plan, avec condition 
de non-glissement a la frontiere. Les donnees initiales que nous considerons sont des 

. perturbations d'energie finie d'un tourbillon regulier dont la circulation a I'infini est 

■ petite, mais nous n'imposons aucune autre restriction a leur taille. En utilisant une 
, estimation d'energie logarithmique et des arguments d'interpolation, nous montrons que 

■ la solution converge lorsque t — >■ oo vers un tourbillon d'Oseen autosimilaire. Ce resultat 
' a ete obtenu en collaboration avec Y. Maekawa (Universite de Kobe). 
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1. Introduction 

(n: 

OO ■ We consider the free motion of an incompressible viscous fluid in a two-dimensional exterior 

*^ domain $7 = \ , where K d^M? \s a compact obstacle with a smooth boundary. We do 

not assume that K is connected, hence we include the case where the fluid moves around a 
finite collection of obstacles. As for the boundary conditions, we suppose that the velocity 
psj \ of the fluid vanishes on dVl and decays to zero at infinity. The evolution of our system is 

thus governed by the Navier-Stokes equations 

{dtu + {u ■ V)n = Au — Vp , div u = , for x e ft , t > , 
u{x,t) = , ioixGdn, t>0, (1.1) 

u(x,0) = uo(x) , for X G , 

where u{x, t) G and p{x, t) G M denote, respectively, the velocity and the pressure of 
the fluid at a space-time point {x,t) G fl x R_|_. As can be seen from the first equation 
in (jl.ip . we assume that the fluid density is constant and that the kinematic viscosity is 
equal to 1. Since p.ip includes no forcing, the motion of the fluid originates entirely from 
the initial data uq : — t- R^, which we assume to be divergence- free and tangent to the 
boundary on 9(7. 

The behavior of the solutions of (jl.ip depends in a crucial way on the decay rate of the 
velocity field u{x,t) as |x| — )• oo. If the initial velocity uq belongs to the energy space 



Ll{n) = {nG L2(J7)2 



div u = in Q , u ■ n = on dfl 
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where n denotes the interior unit normal on dO,, we have the following classical result : 
Theorem 1. For all initial data uq G L'^{0,), Eq. (jl.ip has a unique global solution 

u £ c^{[o,oo),Ll{n)) n c\{o,oo),Ll{n)) n c^i{o,oo),H^{n)^ n H^{nf) , 

which satisfies for all t > the energy equality : 



If 1 

"^(•'Olli2(f7) + / l|V'u(-,s)|||2(f^)ds = -||'Uo||i2(f^) . (1.2) 

~ i/ 

Global well-posedness for the Navier-Stokes equations was first established by Leray 
in the particular case where = M^. When C is bounded, the first results 
also go back to Leray |2S] , but global existence of large solutions was shown only later by 
Ladyzhenskaia [21] , see also [271 US M ■ To prove Theorem [H one can use a regularization 
or a discretization procedure to construct global weak solutions of (jl.ip which satisfy the 
energy inequality, and then prove that these solutions are unique and have the desired 
regularity. Alternatively, one can construct local mild solutions by transforming (jl.ip into 
an integral equation and solving it by a fixed point argument, and then use the energy 
equality (jl.2p to show that all solutions can be extended to the whole time interval [0, oo). 
Although most of the literature is devoted to the situation where is either a bounded 
domain or the whole plane M'^, the case of an exterior domain can be treated without 
essential modifications, see e.g. [22] 



It follows from ()1.2p that the kinetic energy E{t) = ^||'u(-, t)||^2(Q) is nonincreasing in 
time, and a result of Masuda [28] shows that E{t) converges to zero as t — )• oo. Moreover, 
under additional assumptions on the initial data, it is possible to specify a decay rate in 
time. For instance, if uq S L'^{Q,) n L'^{Q,)'^ for some q S (1,2), the solution of (jl.ip lies in 
the same space for all t > and 

II^^(->*)IIl2(c) = as t^oo, (1.3) 

see [HI^IKI]- It is interesting to notice that (|1.3p fails in the limiting case q = 1- Indeed, if 
uq £ L^(r2)nL-^(r2)^, then in general the velocity field u{x,t) decays like as — )• oo, 
so that u{-,t) ^ L^{n)'^ for t > 0. 

As an aside, we mention that this loss of spatial decay is related to the net force F 
exerted by the fluid on the obstacle K. To see this, we flrst observe that any velocity field 
u G L^(0) n L^(r2)^ satisfies f^udx = 0. Indeed, if u is smooth and compactly supported, 
then using Gauss' theorem and the fact that u • n = on we find 

/ Uj dx = / {u ■ V)xj dx = / div(tiXj)dx = , for j = 1,2 . 
Jn Jn Jn 

The general easily case follows by a density argument [20]. Now, if n G C^{[0,T], L'^{nf) 
is a solution of the Navier-Stokes equation (jl.ip . then 

= [ udx = [ (Au -Vp- {u- V)n) dx = - f (Tn) da = -F , 
di Jn Jn^ ' Jan 

because Auj — d-ip = djTij where Tjj = diUj + djUi — p5ij is the stress tensor (we recall 

that all physical parameters have been normalized to 1). The formal calculation above 

can be made rigorous ^20j and shows that, no matter how localized the initial data may 

be, the velocity field u{-,t) does not stay integrable for positive times, unless the net force 

F vanishes identically. Of course this is not the case in general, but in highly symmetric 

situations it is possible to construct solutions of (jl.ip for which F = 0, and which decay 

faster as t — oo than what is indicated in (jl.3p . see [15| [T6]. 

Much less is known about the solutions of (jl.ip if the initial data uq are not square 

integrable. Although the physical relevance of infinite-energy solutions can be questioned, 

we believe that such solutions naturally occur when studying the dynamics of (jl.ip in a 
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two-dimensional exterior domain fi. One way to realize that is to consider the relation 
between the velocity field u and the associated vorticity uj = diU2 — d2Ui- Given p £ [1, 2), 
we denote 

WQ^^{n) = |n G L'^inf Vu G LP{n)^ , dwu = Omn, li = on dQ^ . 

In other words M^qctC^) is the completion with respect to the norm u i— )• ||Vu||lp of the 
space of all smooth, divergence-free vector fields with compact support in fi, see [lUj . We 
then have the following result : 

Lemma 2. If u € Wq'^{Q,) for some p G [1, 2), and if ui = diU2 — d2Ui, then 

= / i^^^(y)dy , (1.4) 

for almost every a; G il. Here = {—X2,xi) and |xp = x\ + X2 if x = (xi,X2) G M^. 

The proof of Lemma [2] is very simple : if n : — )• denotes the extension of u 
by zero outside VL, then u G L2p/(2-p)(m2)2^ ^ LP{M?f, and divu = 0. Moreover 
u) = diU2 — d2Ui is the extension of u by zero outside Q. Thus u can be expressed in 
terms of (D using the classical Biot-Savart law in M?, and restricting that relation to Q we 
obtain (jl.4p . We emphasize that the representation (|1.4p is only valid if to is the curl of 
a divergence- free velocity field u which vanishes on dQ. In contrast, if w : — )• M is an 
arbitrary smooth function with compact support, the velocity field defined by (jl.4p does 
not even satisfy u ■ n = on ! 

We now assume that the vorticity distribution u is sufficiently localized so that u) G 
L^{Q,), and we define the total circulation 

a = uj{x)dx = lim (p ui da^i + U2 dx2 , 
Jn R-^°°J\x\=R 

where the second equality follows from Green's theorem, since uj = diU2 — d2Ui and u 
vanishes on dVl. Using the vorticity formulation of the Navier-Stokes equation (|l.ip . it is 
not difficult to verify that the total circulation is a conserved quantity. But it follows from 
(fT^ that 

u{x) ~ — - — 7K , as Ixl — )• 00 , (1.5) 
27r 

hence u ^ L^iVLY as soon as a 7^ 0. This shows that finite-energy solutions of (jl.ip 
necessarily have zero total circulation. On the other hand, in many important examples of 
two-dimensional fiows such as vortex patches, vortex sheets, or point vortices, the vorticity 
distribution typically has a constant sign, hence the total circulation is necessarily nonzero. 
In our opinion, it is thus important to enlarge the class of admissible solutions of (jl.ip . 
so as to allow for velocity fields which decay like \x\~'^ as |x| — )• 00. 

A possible framework for the study of infinite-energy solutions of the Navier-Stokes 
equation (jl.ip is the weak energy space 



Ll'°°{n) = {uG l2'°°(17)^ 



div u = in i7 , u • n = on 3^2 



where L2'°°(J7) is the weak space on 17, see [3]. We recall that 

/ \ 1/2 

||M||i2,oom) w sup Al meas{x G O I |u(x)| > A} I , (1-6) 
A>o ^ ^ 

in the sense that the norm ||n||^2,oo is equivalent to the quantity in the right-hand side of 

(USD. Clearly LKQ) ^ L^'°°(17), but the weak energy space is large enough to include 

velocity fields which decay slowly at infinity, as in (jl.Sp . Concerning the solvability of (jl.ip 

in Lo-'°°(ri), the following general result was obtained by Kozono and Yamazaki : 



3 



Th. Gallay 



Theorem 3. [23j There exists e > such that, for all initial data uq G La'°°{^) satisfying 

limsup A( measjx G 0| |uo(3;)| > A} j < e , (1.7) 

A-i>+oo ^ ' 

Eq. (jl.ip has a unique global solution such that, for all T > 0, 

sup ||'u(-,t)||i2,oom) + sup t^/'^\\u{-,t)\\L4r^) < OO , 
0<t<T 0<t<T 

and such that u{-,t) uq as t ^ in the weak-* topology of La'°°{i^)- 

Theorem [3] shows that the Cauchy problem for the Navier-Stokes equations (jl.ip is 
globally well-posed in the weak energy space Lo-'°°{Q), provided that the local singularity 
of the initial data uq is sufficiently small, in the sense of ()1.7p . To illustrate the meaning 
of this smallness condition, we consider the simple situation where the initial flow is just 
a point vortex of circulation a G M located at xq G The initial vorticity is thus given 
by uJo{x) = a6{x — xq), and using the classical Biot-Savart law in the exterior domain Cl 
(see e.g. |17] ) it is easy to verify that the corresponding velocity field uq lies in La''^{fi), 
is smooth in J7 \ {xq}, and satisfies 



a {x — xq)- 
27r \x — xo\ 



Uo{X) W — "i [2~ X — )• Xo 



so that (|1.7p is fulfilled if and only if |a| < -v/ivre. This example shows that, if the initial 
vorticity uq is a finite measure, condition (II. 7p implies a restriction on the size of the atomic 
part of ujQ. Such a restriction also arises in the analysis of the two-dimensional vorticity 
equation in the whole space M^, see [13], but in that particular case the uniqueness of the 
solution can be established when the initial vorticity is an arbitrary finite measure |1H [2]. 

Although Theorem O provides the existence of a large class of infinite-energy solutions, 
very little is known about the asymptotic behavior of these solutions as t — >• oo. In fact, we 
do not even know whether they stay bounded in the weak energy space because 
we are lacking a priori estimates. Indeed, if uq ^ L'^{^) the energy equality ()1.2p does not 
make sense, and because of the no-slip boundary condition on dfl it is quite difficult to 
obtain estimates on the vorticity distribution if 7^ M^. In the rest of this paper, however, 
we consider a particular class of infinite-energy solutions of the Navier-Stokes equations 
()l.ip . for which the asymptotic behavior in time can be accurately described. 



2. Main Results 



In the particular case where = M^, the Navier-Stokes equations (jl.ip have a family of 
self-similar solutions of the form u{x,t) = aQ{x,t), p{x,t) = a'^Il{x,t), where a G M is 
a free parameter (the total circulation) and 

= , Vn(x,t) = ^|e(x,t)|2 . (2.1) 

2tt |x|^ V / |x|^ 

These solutions are usually called the Lamb-Oseen vortices. If u{x,t) = aQ{x,t), the 
corresponding vorticity distribution is uj{x,t) = aE{x,t), where 

E{x,t) = die2{x,t) - d2eiix,t) = ] . (2.2) 

47r(l + 1) 

Note that H(a;, t) > and /jj2 '^{x, t) dx = 1 for all t > 0. Oseen vortices play an important 
role in the dynamics of the Navier-Stokes equations in M^. In particular, we have the 
following result : 



4 



THE NAVIER-STOKES EQUATION IN AN EXTERIOR DOMAIN 



Theorem 4. [13]. For all initial data uq G L'^'°°{E?) such that the vorticity distribution 
ujQ is integrable, the solution of the Navier-Stokes equation in M? satisfies 

/ \uj{x,t) — aE{x,t) \ dx > , where a = ujodx . 

In other words, Oseen vortices describe the leading order asymptotics of all solutions of 
the Navier-Stokes equations in with integrable vorticity distribution and nonzero total 
circulation. 

In the case of an exterior domain = M? \ K, approximate Oseen vortices can be 
constructed in the following way. Let x : — ^ [0, 1] be a smooth, radially symmetric cut- 
off function such that x is nondecreasing along rays, x = on a neighborhood of K, and 
xix) = 1 when |x| is sufficiently large. The truncated Oseen vortex with unit circulation 
is defined as follows : 

u^{x, t) = x{x)e{x, = ^ ^ (l - e-^)x{x) . (2.3) 

For any t > 0, it is clear that u^{-,t) is a smooth divergence-free vector field which vanishes 
in a neighborhood of K. The corresponding vorticity distribution oj^ = diU2 — d2u\ has 
the explicit expression 

uj^{x,t) = x(a;)H(x,t) + ^^(l-e"W))x-Vx(a;) , (2.4) 

where H(x,t) is defined in ()2.2p . In particular ijj^{x^t) > and /jg2 (^-^(x, t) dx = 1 for all 
t > 0. Moreover, a direct calculation shows that 

{u^ . V)u^ = -Vjn'^l^ + (n^)^a;X = -^\u^\^ , (2.5) 
2 |xp 

hence there exists a radially symmetric function p^{x,t) such that — Vp-^ = {u^ ■ 'V)u^. 
Now, given a S M, we consider solutions of (jl.ip of the particular form 

u{x,t) = au^{x,t) + v{x,t) , p{x,t) = a'^p^{x,t) + q{x,t) , (2.6) 

where u^{x, t) is the truncated Oseen vortex defined in ()2.3p . and v{x,t) is a finite- energy 
perturbation. In this situation, we expect that v{-,t) converges to zero in energy norm as 
t — )• oo, so that the long-time behavior of u{-,t) is described, to leading order, by the Oseen 
vortex Q0(-,t). Our main result, which was obtained in collaboration with Y. Maekawa, 
shows that this is indeed the case, provided the total circulation a is sufficiently small. 

Theorem 5. [12j Fix q € (1, 2), and let ^jl = l/q—\/2. There exists a constant e = e{q) > 
such that, for any smooth exterior domain Q C M? and for all initial data of the form 
uq = au^{-,0) + vq with \a\ < e and vq € L'^{^) H L^(ri)^, the solution of the Navier- 
Stokes equations (jl.ip satisfies 

\\u{;t) - aQ{;t)\\L^(n) +t'^'\\^u{;t) - aVe{;t)\\L^(n) = , (2.7) 

as t oo. 

To understand the scope and the limitations of this statement, a few comments are in 
order. First of all, Theorem [5] is a global stability result for Oseen vortices with small circu- 
lation at infinity, because we do not impose any restriction on the size of the perturbation 
Vq G L'^i^) n L'^(rj)^. In the particular case where a = 0, there is no vortex at all and we 
just recover the asymptotics (jl.3p with 0(t~^) instead of o(t~^) in the right-hand side. 
Also, in the simple situation where $7 = M^, our result is comparable to that of Carpio [5], 
although the proof is very different. 

The main limitation of Theorem [5] is of course the restriction on the size of the cir- 
culation a, which we believe is purely technical. In this respect, the fact that e{q) can 



5 



Th. Gallay 



be taken independent of the domain is quite significant, because we know that there 
is no restriction on the circulation in the particular case where = M^, see Theorem [H 
Obviously e{q) is a decreasing function of q, and the proof shows that e{q) = 0(\/2 — q) 
as (7 —7- 2. Thus the limiting case q = 2 is not included, which means that we are not able 
to control arbitrary finite-energy perturbations of the Oseen vortex (see however [18] for 
a partial result in that direction). On the other hand the limit of e{q) as (7 — t- 1 can be 
estimated and is found to be approximately e,, = 5.306, see [12J. 

We also mention that the decomposition uq = au^{-,0) + vq of the initial data is 
automatically satisfied if we assume that the initial vorticity is sufficiently localized. More 
precisely, we have the following auxiliary result, which follows quite easily from Lemma [2l 

Proposition 6. [12J Given q G (1,2), assume that uq S Vi^Q'^(rj) for some p G [1,2) and 
that the associated vorticity ujq = curl uq satisfies 



for some m > 2/q. If we denote a = ji^ujQ{x) 6.x, then uq = au^{-,0) + vq for some 
vq £ Lai^) nL'^(il)^. In particular, if \a\ < e, the conclusion of Theorem\^ holds. 

In view of Proposition [6l it would be nice to extend the conclusion Theorem [5] so as 
to include a convergence result for the vorticity distribution in the critical space L^{Q,). 
This is not immediately obvious, because the classical L^-L'^ estimates for the Stokes 
semigroup in an exterior domain do not include the limiting case p = 1, see [H^. However, 
combining Theorem [S] with a relatively standard estimate, which shows that the norm 
of the vorticity cannot leak to infinity, we obtain the following result which is the main 
original contribution of the present paper : 

Proposition 7. Under the assumptions of Theorem\^ if we suppose in addition that (|2.8|) 
holds for some m > 2/q, then the vorticity uj = curln satisfies 



where E{x,t) is defined in (|2.2p . 

In the rest of this paper, we give a simplified proof of Theorem [5l which does not yield 
the optimal conclusion. In particular, we shall find a suboptimal convergence rate in (|2.7p . 
and our limitation on the size of the circulation will (a priori) depend on the domain il. 
We refer the reader to for a complete proof, including all details. In the last section, we 
briefly show how Proposition [7] follows from Theorem [5l using some additional information 
on the vorticity near infinity. 

3. Energy estimates 

Given a G M we consider solutions of (jl.ip of the form ()2.6p . The perturbation v{x,t) 
vanishes on the boundary dil. and satisfies the equation 

dtv + a{u^ ■V)v + a{v ■V)u^ + {v ■V)v = Av + aR>^ -Vq , divv = , (3.1) 

where is the remainder term given by (|6.5p below. It is not difficult to verify that the 
Cauchy problem for equation (jS.ip is globally well-posed in the energy space -L^(il). The 
goal of this section is to control the long-time evolution of the perturbation v{t) = v{-,t) 
using energy estimates. 

First of all, we multiply both sides of (j3.ip by v and integrate by parts over il. Taking 
into account the no-slip boundary condition, we find 




(2.8) 




as t — 00 



(2.9) 




^;(t)||i. + ||V^;(t)||i. 



a{v{t),R^{t)) - a{v{t), {v{t) • V)u^{t)) 



(3.2) 
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where (•, •) denotes the usual scalar product in L^(il), so that || • ||^2 = (• , ■^^'^ . To estimate 
the right-hand side of (j3.2|) . we use the results of Section [U] below. First, in view of (|6.7p . 
we have 



lal ,,,, /x,,9 Kq 

J^IIV.WII., <-||V„(t)||i, + ^|JL_ 



\a{v(t).n^m < ^||V..(t)||i, < -||Vi. 
Moreover, applying (|6.2p with p = 00, we see that 



\{v{t),{v{t).V)uHt))\ < ^ll^Wlli. 



b. 

1 + t 

We thus obtain the energy inequality 

^^imwh + iiv^^wiii. < ^ibwiii. + , t>o. 

Using Gronwall's lemma, we deduce that 

IK<)lli..£l|v.(„),li,d, . (i±l)-'"'(,K,,|,|,.^i^) , ,3.3, 

for t > tQ > 0. This simple estimate shows that the energy of the perturbation v{x,t) 
grows at most polynomially in time as i — )• 00. Such a conclusion is rather pessimistic, 
however, because by a relatively simple modification of the previous argument it is possible 
to establish a logarithmic bound, which is clearly superior for large times. 

Proposition 8. There exists a constant Ki > such that, for any circulation a € M and 
any vq G L'^{Q), the solution of ()3.ip with initial data vq satisfies 

Mml2+ l^\\Vv{s)\\l,ds < 2^^"^Ki(\\vo\\l2+a^log{l + t)) , (3.4) 

for all t >0, where c = 26oo > 0. 



Proof. If t < 1 then (|3.4p follows from p.3p with to = 0, hence we can assume that t > 1. 
Given any r > 0, we denote 

v{x,t) = u{x,t) — au^{x,t + t) = v{x,t) + a(^u^{x,t) — u^{x,t + t)^ , (3.5) 

for all X €z Q and all t>0. Then v satisfies (|3.ip where u^{x,t) and R^{x,t) are replaced 
by u^{x,t + T) and i?^(x, t + r), respectively. Proceeding exactly as above, we thus obtain 
the energy estimate 

\\vmh+ l\\yHs)\\hds < (i±l±I)''"'(||i)(0)||i.+Ca2) , t>0. (3.6) 

Now, we fix t > 1 and choose t = t — 1. From (j6.3p . (j3.5p . we have 

\\v{t)\\l2 < 2\\v{t)\\l2+2a^\\u^{t)-u^{2t-l)\\l2 < 2\\i{t)\\l2 + 2Kia^ log2 . 
Similarly, using (j6.4p . we find 

/ \Nvis)\\l2ds < 2 [ \\Vv{s)\\l2ds + 2a^ [ \\V (s) - Vu^{s + t - ds 
Jo Jo Jo 

<2 f \\Vv{s)\\l2ds + 2K2a^ log . 
Jo ^ 

Thus it follows from (j3.6p (with t = t — 1) that 

\\v{t)\\l2+ l\\Vv{s)\\l2ds < 2=H+i(||^(0)||i.+Ca2)+Ka2log(l + t), (3.7) 

for all t > 0, where k = 2max(K;i, K2). Finally, we have by ()6.3p 

Ih < 2\\vo\\l2+2a'^\\u^{0)-u^{t-l)\\l2 < 2||^;o|li2 + 2Kia2 bgi , 
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hence (|c{.4|) easily follows from (|3.7|) . □ 

Remark 9. The logarithmic energy estimate (j3.4p is the main new ingredient in the proof 
of Theorem To a certain extent, we use it as a substitute for the energy equahty (|1.2p . 
which does not make sense for the solutions we consider. As is clear from the proof, the 
logarithmic energy estimate relies on the fact that Oseen's vortex (j2.ip has "nearly finite 
energy", in the sense that the integral defining ||0(-,t)||^2 diverges only logarithmically 
at infinity. 



4. Fractional interpolation 

Let P be the Leray-Hopf projection in and A = —PA be the Stokes operator, see e.g. 
[U]. We recall that A is self-adjoint and nonnegative in L^(ri), so that the fractional power 
A'' can be defined for all /3 > 0. The following result shows that the range of A^^ contains 
the (dense) subspace LK^I) D L''{Q)'^. 

Lemma 10. H [21] Lei g G (1, 2) and = ^ - \. For all v G Ll{Q) n ^^(O)^, there 
exists a unique w € D{A'^) C L'^{^) such that v = A^^w. Moreover, there exists a constant 
C > (independent of v) such that ||if||L2(Q) < C'll'y||L9(r2)- 

Remark 11. If v,w are as in the above statement, we denote w = A^^''v. Roughly 
speaking, the proof of Lemma [10] argues as follows. By classical Sobolev embedding, we 
know that the domain of A^ is contained in L'^{^) H L'^ (f^)^, where J7 = ^~/^ = 1~|) 
and by duality we deduce that the range of A~^^ contains L'^{^) H L'^iyt)'^. 

We go back to the study of the perturbation equation (|3.ip . which can be written in 
the equivalent form 

dtv + Av + aP[{v^ ■ V)v + {v ■ V)u^^ + P{v ■ V)v = aR^ . (4.1) 

So far we only considered solutions in the energy space -L^(r2), but now we assume in 
addition that vq G L''(J7)^, for some fixed q G (1,2), and we denote /i = | — ^ G (0, ^). 
Then it is not difficult to verify that the solution v{t) of (fO]) lies in Ll{n) n Li{nf for 
all t > 0. In particular, invoking Lemma [TOl we can define w{t) = A~^v{t) for all t > 0. 
This quantity solves the modified equation 

dtw + Aw + aFf,{u^,v) + aF^{v, u^) + Ff,{v, v) = aA-^'R^ , (4.2) 

where F^{u,v) is the bilinear term formally defined by 

F^{u,v) = A-^'P{u-V)v . (4.3) 

We refer to |21l Section 2] for a rigorous definition and a list of properties of the bilinear 
map Ff^ . Our goal here is to establish the following estimate : 

Proposition 12. There exists > and, for all a G M, there exist positive constants 
K2{a) and k{a) such that, if v is any solution of (|4.ip with initial data vq G L'^{^) n 
L^(r2)^, the function w{t) = A~^^v{t) satisfies 

\\wit)\\l,+l\\Vwis)\\l,ds < (l + 0"''=(")exp(if2(a)||^o|li2+i^3)(lbo||L+«'), (4.4) 

for all t > 0. Moreover K2{a) and k{a) are 0(1) as a — )■ 0. 

Proof. Taking the scalar product of both sides of (j4.2p by w, we obtain 

~\\wmh + \\A'/Mmh+a{F,{u^t),v{t))Mt))+»{FM 

+ {F^{v{t),v{t)),w{t)) = a{A-^R^{t),w{t)) . (4.5) 
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It is well known that \\A^I'^w\\l2 = \\Vw\\l2 for all w G D(^V2) = Ll{yi) n HliVtf . To 
bound the other terms, we observe that 

\{F^,{u^,v),w)\ = \{{u^ ■V)v,A-^w)\ = \{{u^ ■V)A-^'w,v)\ 

< ||n^||Loo||A^~''?i;||^2||i;||i2 = ||u^||Loo||^3-/^tt;||^2||74'^w||2,2 

< lln-^llioo II A"^/^zi;||j^2 ||t(;||i2 , 

where in the last inequality we used the interpolation inequality for fractional powers of A. 
The same argument shows that \{F^{y^u^),w)\ < ||u-*^||loo || A-'-/^u;||2,2 ||'u;||j;^2 . In a similar 
way, 

\{F^{v,v),w)\ = \{{v-V)v,A->'w)\ = \{{vV)A->'w,v)\ 

< ||?;|||4p^-''?i;||^2 < C||V7;||i2||t;||i2p5-/^7i;||^2 

< C||Vw||i2p^/2^||i2||u;||i2 . 

Finally, since \{A-~^'R^,w)\ = \{R^ , A">'w)\ < /t3(l + t)-i||yl5-/^u;||^2 by ([621), we can use 
interpolation and Young's inequality to obtain 



for some 71,72 > 1 satisfying 72 + 2/^71 = 2. Thus ()4.5p implies 

^||«;||i2 + ||Vu;||i2 < -||Vu;||i2 + C7||Vu;||i2||«;||i2(|a|||n>^||L^ + ||Vt;||i2) 

1,,^ „2 2||u;||22 2Ca^ 

+ 2ll^-lli^ + (TT#r+(rTI)^ (4.6) 

< CJ^||i2 (a^lln^^lli. + ||V.||i2 + + 



for some positive constants Ci , C2 . 

Now, using ()6.ip with p = 00 and the logarithmic energy estimate ()3.4p . we obtain 

ds 



C, [\a^uHs)\\l^ + ||Vt;(.)||i2 + 7^^) 



(l + s) 

< a^k{a)\og{l + t) + K2{a)\\vo\\l2+C3 , 

where ^2(0) = 2^l"'(7ii^i, k{a) = Cia^ + i^2(a), and C3 = Ci(7i - 1)-^ Applying 
Gronwall's lemma to (I4.6I1. we thus find 



Mmh + ||V«;(.)||i2ds < (l + tr'*=(")exp(K2(a)lbo|li2+C73)(|ko||i2+C4a2) , 

where C4 = 02(72 — 1)"^, and (j4.4p follows since ||wo||l2 < C||7;o||l9 by Lemma [TUl □ 

Corollary 13. Under the assumptions of Proposition [7^ there exists a positive constant 
K4 depending on |a| and Ht'ollL^nL'J such that, for any t > 2, there exists a time to E [t/2, t] 
for which 

\Hto)\\h < i^4(l + to)"''(")''^ . (4.7) 
Proof. Fix t > 2. In view of (|4.4p . there exists a time to £ [V^)^] such that 

||Vu;(to)||i2 < y r ||Vu;(s)||i2d5 < ^ Ko(l + 1)"'^(") < 2"'^(")+2i^o(l + to)"''^")"' , 

I Jt/2 t 
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where i^o = exp(i^2||t'o|li2 +i^3)(||^'o|li. Moreover, \\w{to)\\l, < Ko(l + to)"'^(") by 
()4.4p . Thus, we obtain ()4.7p using the interpolation inequaUty ||t;(to)||L2 = II^'^^(^o)||l2 < 
\\Vwito)\\f,\\wito)\\]^'^. □ 

We are now able to conclude our sketch of the proof of Theorem [5j Given a G M and 
vq £ L'^i^) n L'?(r2)^, let v{x,t) be the solution of the perturbation equation (|4.1|) . For 
any t > 2, we choose Iq G [t/2,t\ as in Corollary [T31 and we apply estimate (|3.3|) . We thus 
obtain 

\\v{t)\\l, < 22^-H(||t;(to)||i2 + ^^) < C^ll+tr'^^")-^'^ , (4.8) 

where Ca > is 0(1) as a — )■ 0. Now, if |a| is small enough to that a^ki^a) < 2fi, the 
right-hand side of ()4.8p converges to zero (at a suboptimal rate) as t — )■ oo. In particular, 
the perturbation v{-,t) becomes very small in energy norm for large times. In that regime, 
the perturbation equation (|4.ip can be solved by a global fixed point argument, which 
allows to show that 

H;t)h2^n)+t'^^\Nvi-MLHn) = Oit'^ , (4.9) 

as t ^ oo, see [121 Section 3]. Finally (|2.7p follows from (j4.9p . because v{x,t) = u{x,t) — 
au^ix,t) and ||u^(-,t) - e(-,t)||i2 + ||Vn>^(-,t) - VG(-,t)||i2 < C(l + t)"^ for aU t > 0. □ 



5. Convergence of the vorticity 

This final section is devoted to the proof of Proposition [71 We first show that, under the 
assumptions of Theorem [5l one can control the norm of the vorticity sufficiently far 
away from the obstacle K. 

Proposition 14. Under the assumptions of Proposition^ the vorticity uj = curln satisfies 

\uj{x,t)\dx = Oit-^") , as t^ CO . (5.1) 



i 



|x|>tl/2logi 



Proof. Since by ()2.8p the initial vorticity is assumed to be square integrable, the solution 
u{x,t) of (jl.ip given by Theorem [5l satisfies 

||n(.,t)-aG(-,t)||i2(o) + (l+t)^/2||Vn(.,t)-aVe(-,t)||i2(f,) < , (5.2) 

for all t > 0, where Co > depends only on the initial data. The associated vorticity 
oj = curl n is a solution of the advection-diffusion equation 

dtuj + u- Vuj = Auj , X , t > , (5.3) 

but the no-slip boundary condition becomes very complicated when expressed in terms 
of oj. It is thus difficult to use (j5.3p to obtain estimates in the whole domain fi, and in 
particular near the boundary 90. Here, however, our goal is to bound uj near infinity, so 
we can avoid that problem using localized energy estimates and invoking (|5.2p to control 
the fiux terms in the regions where the localization function is not constant. 
Given T > 4 and i? > 1, we define the cut-off function 

where r{t) = 2^3/2^1/2 iog(t/2), and : [0, oo) [0, 1] is a smooth, nondecreasing function 
satisfying (/)(r) = for r < 1 and <j){r) = 1 for r > 2. We always assume that T > 4 is large 
enough so that the support of ip{-,t) is contained in 0,, and that i? > 1 is large enough 
(depending on t and T) so that ip{-,t) is not identically zero. Given A > 0, we also denote 
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and we observe that < ^x{lo) < \uj\ and ^'{{lo) > for all w G M. 
Now, using (j5.3p . we obtain by a direct calculation 

^ / ?/;$A(w)dx = / + AV' + (n • V)V')$A(w)dx - / V^A(a;)|Va;|^ dx 
Jn Jn^ ' Jn 



because ipt 1^ and ^^i^) > 0. If we integrate this inequality over t G [0,T], we find 

^(x,r)^>A(a;(2;,T))dx < [ i;{x,0)<^>x{uJo{^))d^+ [ [ Q,^,«(x, t)$A(w(x, t)) dx dt , 

Jn Jo Jn 

where Q^^u = I^V' + " Using Lebesgue's monotone convergence theorem, we can 

take the limit A — t- in both sides, and we arrive at the simpler estimate 

'ip{x,T)\u!{x,T)\dx < / ■iP{x,0)\ujq{x)\ dx + / / Q^^u{x,t)\uj{x,t)\dxdt . (5.5) 

Jn Jo Jn 

Our next task is to take the limit — t- oo in (|5.5|) . Again, we use the monotone 
convergence theorem, except in the last integral where it does not apply. To treat that 
term, we observe that Q^^ui'ii) vanishes identically except in the region Dr U 
where for any /> > we denote Dp = {x G | p < |x| < 2p}. Taking R > sufficiently 
large and using ()5.2p . we easily obtain 

/ Qi,,u\uj\dx < [ (\AiP\ + \u\\ViP\)\uj\dx < ^ , 
Jdr Jdr^ ' R 

for some Ci > independent of t G [0, T]. The contribution of the annulus Dfi is therefore 
negligible for large R, hence taking the limit -R — )> oo in (|5.5p we arrive at 

/ iJi{x,T)\uj{x,T)\dx < / iJi{x,0)\ujo{x)\dx + / / Q^,,„(x, t)|a;(x, t)| dx dt , 
Jn Jn Jo Jn 

where 'ipi{x,t) = (j){\x\/r{t + T)). In particular ^i(x,T) = 1 for |x| > T"^/^logT and 
'tpi{x,0) = for |x| < r(T), hence the last inequality implies 

/ |a;(x,r)|dx < / |ti;o(x)|dx+ / / Q^^^u{x,t)\uj{x,t)\dxdt . (5.6) 

y|x|>TV2iogT J\^\>r{T) Jo Jn 

To conclude the proof of Proposition [T^ it remains to estimate both terms in the 
right-hand side of ()5.6p . First, using ()2.8p and Holder's inequality, we easily find 

1/2 / „ . \ V2 



|a;o(x)|dx < ( / (l+|x|2)-|^o(x)|2dx) ( / ^_^^dx 



1 

x\>r{T) \Jn J \J\A'>i'^) 

(l+|x| 

< C7r(T)-(™-i) < C^T-^" , (5.7) 



for some C2 > independent of T . In the last inequality, we used the hypothesis m > 
Ijq = 1 + 2/2 and the fact that r(T) > CT^/'^. On the other hand, since ^i(x,t) is given 
by ()5.4p with R = 00, there exists C3 > such that 

|vV'i(x,t)| < ^:^^^D', , \AMx,t)\ < , 

where D'^ = = {x G | r{t + T) < \x\ < 2r{t + T)}. It follows that 

J^Q^U^,t)Hx,t)\dx < C3 H^,t)\dx. (5.8) 
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But using (j5.8|) and (|5.2|) we easily find, 

uj{x,t)\dx < / oj{x,t) — a'B{x,t) d2; + |a| / E{x,t)dx 



< measpO ^ " aS(-, Olli^cn) + Ml^i', t)\\LHD',) 

< c/^^!TL+cexpf--(;+^^;) < c 



(l+t)M+l/2 ' 4(l + i);- (l + t)M+l/2 ' 

where in the last inequality we used the fact that r{t + T) > C{t + T)^^^ log(t + T). In a 
similar way, 



\u{x,t)\\uj{x,t)\dx < / \u{x,t)\\uj{x,t) — aE{x,t)\dx + \a\ / t)| t)| dx 

, \ C 

< \\u\\LHD',)[\\^i-,t)-aE{-,t)\\L2^n) + M-,t)\\L^D',)) < ' 



Inserting these estimates in the right-hand side of ()5.8p . we obtain 



T 




C C4 
ufa^i t)| dx dt < / — — ^ — TTTTT dt < — — , (5. 



Jn 



for some C4 > independent of T. Thus, if we combine ()5.6p . ()5.7p . and ()5.9p . we conclude 
that 

|a;(x,r)|dx < % , 

|x|>ri/2iogT 

for all sufficiently large T > 0, which is the desired result. □ 

It is now easy to conclude the proof of Proposition [71 For t > sufficiently large, we 
denote fit = {x £ fl\\x\ < t^/^ logt} and we decompose 

uj{x,t) — aE{x,t)\ dx < / \u!{x,t) — aE{x,t)\ dx + / ( |a;(x, t)| + |a|H(x, t) j dx . 
n Jvit Jn\nt ^ ' 

The last integral in the right-hand side is controlled using Proposition [T3] and the ex- 
plicit expression ()2.2p of H(x,t). To estimate the first integral, we use ()5.2p and Holder's 
inequality : 

|w(x,t) -aH(x,t)|dx < V^rt^/^logt ||a;(-,t) -aH(-,t)||i2(m < C . 
Summarizing, we find 

\uj(x,t) — aE(x,t)\ dx < C- — , 
for all sufficiently large t > 0. This concludes the proof. □ 

6. Appendix : estimates for truncated Oseen vortices 

In this appendix, we collect a few estimates on the truncated Oseen vortices (|2.3p which 
are used throughout the paper. We first list a few bounds which follow from (|2.3p and 
(|2.4p by rather straightforward calculations, see |12l Lemma 2.1]. 



Lemma 15. 

1. For any p £ (2, 00], there exists a constant Up > such that 



k^(-,i)llLp(R2) < ^ , t>0. (6.1) 

(l + t)2 p 
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2. For any p G (l,oo], there exists a constant bp > such that 



|V'U^(-,t)||LP(M2) 



< 



1 1 ' 
[l+t)^-p 



t > 



3. There exists a constant ki > such that, for all t,s > 0, 



\u'^{-,t) — u^{-, s) 



Il2 



< Ki 



log 



1 + t 



l + S 



4. There exists a constant k;2 > such that, for all t,s >0, 



\Vu^(-,t) 



Vn'^(-,s)||^2(K2) < t2 



1 



1+t l+s 



(6.2) 



(6.3) 



(6.4) 



Since the truncated Oseen vortex is not a solution of the Navier-Stokes equation, we 
also need a control on the remainder term = Au^ — dtu^ = (Ax)0 + 2(Vx • V)0, 
which has the explicit expression 



i?^(x,t) = @{x, t)Ax{x) + 2 ^ " ^^^^^ ( x^E{x, t) - Q{x, t) 

\x\^ \ 

Lemma 16. There exists a constant > such that, for any p G [l,oo] 



I^^(-,*)IIlp(IR2) 



< 



1 + t 



t > . 



Moreover, for any vector field u G H^^^{M?)'^, we have 



R^{x, t) ■ u{x) dx 



< 



^3 
l + t 



LHD) 



t > 



(6.5) 
(6.6) 
(6.7) 



where D dVl is a compact annulus containing the support ofVx- 

Proof It is clear from ([62]) that \R^{x,t)\ < C{1 + t)-^lD{x) for ah a; G g^^^ ^ > 
and ()6.6p follows immediately. Moreover, we have R^{x,t) = x^Q^{x,t) for some radially 
symmetric scalar function Q{x,t), hence R^{-,t) has zero mean over the annulus D. If 
u G -ff[Qj,(M^)^ and if we denote by u the average of u over D, the Poincare-Wirtinger 
inequality implies 

R^{x, t) ■ u{x) dx 



D 



R^{x,t) ■ {u{x) - u)dx < C||i?^(-,t)||i2(R2)||Vn||i2(£,) , 



and using ()6.6p with p = 2 we obtain ([6 



□ 
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